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Abstract 

The relevance of the electronic swap in the stopping process of proton by hydrogen is investigated. 
To this end, the Classical Trajectory Monte-Carlo method is used to calculate the k-stopping cross- 
section, i.e. the stopping cross-section given the occurrence of k-swaps during the collision. It is 
found that electron swaps can be used to label electron trajectories, as it seems to describe fairly 
well the extent of the electron-ion interaction during the collision. Depending on the ion energy, 
the number of swaps entering the stopping cross section may vary. In the keV range the number of 
swaps can be as large as five or more, whereas, at larger energies, only three or less electron swaps 
may take place in collisions of relevance to stopping. 
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I. INTRODUCTION 



The so-called stopping power, or linear energy deposition for swift ions penetrating a 
material, i.e. 5E/<5x, is a relevant quantity in many areas of engineering and physics 
It is normally defined as the mean energy-loss (5E) a projectile undergoes when passing 
through a £x-thick layer of material. From an atomistic point of view the stopping power is 
obtained using the equation: 

where N is the atomic density of the material, and S is the stopping cross-section, defined 

as 



S = J dTa{T) . (2) 

where cr(T) is the differential scattering cross section as a function of the energy loss T. 
The stopping cross section has been extensively studied both experimentally and theo- 

rniflnnna 

retically PJJ2I, y, m, 15|, la, 131 • F rom the theoretical side, following the seminal works of Bohr 
0, Bethe J and Lindhard |4j, a considerable progress has been achieved over the years. 
Furthermore, the increase of computer capability observed in the last decades brought on 
the possibility of numerically simulating the scattering processes, and triggered the interest 
of many researchers in re-studying the stopping at a level of detail not performed before. 
Computer simulations of the stopping are based on both quantum j^, [t| and classical 



mechanics 
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Hj. On the classical side, the so-called Classical Trajectory Monte-Carlo 
method (CTMC) assumes that the nuclei and the electron are all classical particles. The 
simulation, therefore, first picks the initial state of the electron and the nuclei and, secondly 
solves the Newton's equation of motion for these three particles during the passage of the 
bombarding ion. The use of classical mechanics, however, is possibly one of the most serious 
shortcomings of this approach. Except for some cases where using classical mechanics is 
seemingly well justified, the results of these simulations may not be expected to be very 
accurate. Surprisingly enough, however, the results of the stopping cross sections calculated 
using the CTMC method showed that this approach is fairly accurate and that one may rely 
on it much more than expected [ll, l^ . 
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The concept of electronic swap was introduced in Ref.[l3|] in analyzing the charge ex- 
change in ion collisions with Rydberg atoms. Although it will be explained in more detail in 
the following section, the electron swap refers to the passage of the electron from one nucleus 
to the other. Obviously, a swap may take place when the two nuclei are close enough so 
that electron is subjected to the forces of the two colliding nuclei. During the collision, the 
classical electron may realize none, one, two or more swaps. It is evident that the number 
of swaps tells one about the trajectory followed by the electron during the collision. There- 
fore, one may conceivably imagine that the number of swaps must be closely related to the 
stopping process as well. 

The present paper aims at analyzing the correlation between the electron swaps and 
stopping of protons by a classical hydrogen atom. To this end, a previously developed 
CTMC computer code HQ is used, but modified so as to enable it to calculate of 
the stopping cross section for trajectories with the same number of swaps separately. This 
paper is organized as follows: in Section [Til the computer simulations code is briefly described. 
Section [TTTI contains the results and discussions of the present calculations. Finally, summary 
and concluding remarks are offered in Section IIV1 



II. CALCULATIONS. 
A. The CTMC method. 



The classical trajectory Monte-Carlo method (CTMC) assumes that both the electron 
and nuclei are classical, point-like, charged particles. The evolution with time of the position 
and velocity of these particles is obtained by solving the Newton's equation of motion, where 
forces are calculated from the Coulomb potential acting between the electron and nuclei, 
and the initial position and velocity of the electron are randomly chosen using a certain 
statistical initial distribution function. 

An essential part of the CTMC calculation is the one trajectory event or, simply, a 
trajectory. It comprises the preparation of the initial state and the numerical calculation 
of the electron and nuclei trajectories up to a time when the state of the electron can be 
unambiguously determined. That is, the electron has to be exclusively excited, ionized or 
captured by the projectile. 
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A trajectory is completely defined by the initial state of the electron and the impact 
parameter b of the incoming ion. Therefore, any quantity one would like to calculate should 
be obtained by taking average over a certain number of trajectories. A number that often 
depends on the statistical uncertainty one is prepared to accept. For example, the mean 
energy loss undergone by the incoming ion for a given impact parameter b and bombard- 
ing energy e, i.e. AE(b,e), is obtained from calculating the average energy losses over N 
trajectories events as: 

AEJb,e) 

AE(b,e)= ^= 1 N A ) (3) 

where AEj(b, E) is the energy transferred from the projectile to the electron and target 
nucleus in the j-th trajectory event, and which is often obtained from the difference between 
the initial and final energy of the electron. Similarly, one can obtain the stopping cross- 
section by calculating the integral: 

POO 

S(e) = 2n / dbbAE(b,E). (4) 
Jo 

Using the CTMC, however, the integral above is calculated by taking average of 
2-7T b AE(b, E) over a number M of impact parameter b, uniformly sampled between and 
b r . 



J max • 



V*' bj AEj (bj , e) 
S(s) - 2^ W Z ^- 1 3 M 3 3 , (5) 

where both b max and M are chosen from a compromise between accuracy and computer 
time. Sometimes though, b max is obtained from physical arguments. This is the case of high 
bombarding energies where Eqs. (4-5) do not converge unless a finite b max is used. Normally 
bmax is set to the so called adiabatic radius, defined as the impact parameter above which a 
real, quantum electron will no longer absorb energy from the passing charge 



B. Initial distribution. 



Unfortunately enough, the choice of the initial state of the electron is not as simple as one 
could possibly imagine. The initial position and velocity of the electron are sampled from a 
classical phase-space distribution function that should be similar to that of the target atom. 
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However, this is not so simple, since one needs both the position and velocity, and quantum 
mechanics forbids the simultaneous knowledge of these two quantities. The introduction of 
some approximation is therefore unavoidable. 



For example, one can use the so-called microcanonical distribution 
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16| . In this ap- 



proach, the initial state of the electron in the classical orbit is sampled from the distribution: 



p(r,v) 



where S is the Dirac delta function 



87H V 2m r 



(6) 



17|, m is the mass of the electron, e is the elementary 



charge, k e is the so-called Coulomb's constant, and Eq the energy of the ground state in the 
hydrogen atom, i.e. -12.6 eV. This distribution has the advantage of being stable over time 
and that every sampled electron will have the same energy. However, it is limited along the 
radial coordinate and this may cause some underestimation of the charge transfer and the 
stopping cross sections. 



C. Swap number 

As was already mentioned, the electron is subjected to the combined Coulomb interaction 
of the two nuclei, i.e.: 

e 2 e 2 

V e (v) = -k e - r - k e] r , (7) 

|r e -r x | |r e — r 2 | 

where r e is the position of the electron and ri, T2 are those of the ion and target nuclei 
positions, respectively. 

In order to identify under which nucleus the electron is much stronger interacting with, 
the space can be divided into two by a plane that is perpendicular to the straight line 
containing the two nuclei, that passes through the so-called saddle point [18] defined by the 
expression: 

(8) 



|r e -r x | |r e — r 2 | 

In the following, this plane will be referred to as the mid-plane (see figure [1]). One may 
thus introduce the swap number, i.e. Nj, as the number of crossings undergone by the 
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FIG. 1: Path of the electron and position of the mid-plane (crossed stripe) for a v = 1 atomic units 
(a.u.) proton (left) impinging on a Hydrogen atom (right). 

electron through the mid-plane during the j-th trajectory event. Similarly, one can define 
the probability of having fc-swaps after iV trajectories as: 



N 



Pk(b,e) 



where S^j is the Kroenecker delta function 



(9) 



17| . It can be easily verified that 



(10) 



k=0 



And, similarly, the one can introduce the swap number cross-section as: 

POO 

a{e) = 2tt / dbbP k (b,e) . (11) 
Jo 

This cross section is expected to be finite only for k > 0, since the integral above diverges 
for k — 0, because Po(b, e) = 1 with b — > oo. 

Similarly, one may define the &-swap stopping cross-section or, simply, the fc-stopping, as 

POO 

S k (e) = 2n dbbAE k (b,e), (12) 
Jo 

where AE k (b, e) represents the mean energy-loss associated with fc-swap trajectories. 
It must be noticed that the total stopping cross section becomes 




S(e) = ^S k (e). (13) 

k=0 

It must be observed that, contrary to the swap-number cross section, the zero-swap 
stopping cross section can be readily calculated, since AE Q (b, E) does not have any additional 
problem regarding convergence, apart from those pointed out above at the end of subsection 
2a. 

In the following section the results of calculating the swap-number and Sk for protons on 
Hydrogen will be presented and discussed. 

III. RESULTS AND DISCUSSIONS. 

Figure [2] shows the swap number cross-sections Cfc(e) as a function of the bombarding 
energy e and swap numbers ranging from 1 to 7. Notice that the (k = 0)-case does not 
appear since, as was previously commented, a (e) cannot be calculated. In the first place 
one can see that, with the only exemption of k — 2, is a decreasing function the swap 
number. Secondly, ak{e) appears to be a relatively smooth function of the energy for k = 2, 
whereas it strongly depends on e for higher A;-values. Third, it must be noticed that the k = 
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FIG. 3: fc-swap stopping cross-sections [see Eq. (ll2p ] as functions of the bombarding energy. 

2 case remains fairly constant over the entirely range of bombarding energy studied in this 
paper. This implies that the probability for the electron of performing one swing towards 
the projectile and coming back to the target nucleus is nearly the same irrespective of the 
bombarding energy. It is evident that with a decrease of the bombarding energy the number 
of swaps increases. This agrees with the commonly accepted picture that, at low energy, the 
electron performs several turns around the two nuclei during the collision. 

The fact that o<i remains at nearly a constant value with an increase of the ion energy 
is not at all unexpected. As a matter of fact, excepting for an exceedingly small number of 
head-on encounters between the projectile and the electron, the latter is practically at rest 
compared to the projectile velocity and so, the plane used to count the swaps may pass either 
zero or two times over the electron, irrespective of the ion energy. The case of zero-swap 
however is not shown in the figure due to its ill- definition, as was previously mentioned. 

The results of calculating the k-stopping cross-section [see Eq. (ll2|) ] for k = to 5 are 
depicted in figure [31 There one can see that for bombarding energy lower than 100 keV, 
Si dominates, whereas, for larger energies, S2 becomes the largest of them all. It must 
be noticed that both S4 and S5 drop down very rapidly at energies greater than 10 keV, 
becoming negligible small for energies greater than 100 keV. This is due to the fact that 
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FIG. 4: Accumulated fc-swap stopping as functions of the bombarding energy. Sj is obtained from 



swap numbers greater than 3 are fairly unlikely events for energies larger than 10 keV, as 
was already observed in figure [2j 

FigureHJshows the sum of the first k-swap stopping, normalized to the total stopping cross 
section. This plot shows in a much clear manner the contribution of the first /c-stopping to 
the total one. As one can readily see, for energies smaller than 3 keV the number of swaps 
associated with the stopping process is fairly larger than 5. As the bombarding energy 
increases this number decreases, too. However, at large energies the contribution of one- 
and zero-swap collisions decrease and that of the two-swap collisions becomes the dominant 
contribution. 

One may obtain the mean energy-loss in a fc-swap collision by dividing the /c-stopping by 
the /c-swap cross section. The results are shown in the figure El where even and odd swaps 
collision are plotted separately. As one can readily see, there is a clear difference between 
even- and odd-swaps collisions. In the first place, the mean energy-loss corresponding to odd- 
swaps are greater than those of even-swaps. Secondly, the mean energy-losses of odd-swaps 
collisions are monotonously increasing functions of the projectile energy, whereas those of 



Eq. (fT2|) and S is the total stopping cross section, i.e. S = Yl'jLo 
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FIG. 5: Mean energy-loss per collision, denned as the ratio between the /c-stopping and the cor- 
responding A;-swap cross section. Results are grouped by (a) even- and (b) odd-swap collisions, 
respectively. 



the even swaps are bell-shaped with maximum around 100 keV. 

In figure El the total stopping cross-section calculated in this paper is compared with 

1 Q, 



experimental data from Refs. 19 



20 
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23]. As was previously discussed [UJ. 



the present computer simulations reproduce experiments fairly well over a wide range of 
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FIG. 6: Comparison between present calculations (lines) and experimental data (symbols). 



Refs.[19|, 



3, 



21 



22, 



23]. 



bombarding energy. Discrepancies are observed below 10 keV, however, such differences 
may be conceivably attributed to the fact that molecular, not atomic hydrogen targets are 
used in those experiments. With regard to the electron swaps, it must be noticed that 
by adding the first three swap-stopping cross sections, i.e. 0, 1 and 2, one may account, 
reasonably well, for the total stopping. 

IV. SUMMARY AND CONCLUDING REMARKS. 

By using the Classical Trajectory Monte-Carlo method, the correlation between the so- 
called electron swaps and the stopping cross section is investigated. The electron swap, a 
concept introduced in Ref.fisl] in connection with the charge exchange during collision of 
ions with Rydberg atoms, relates to the passage of the classical electron through the ideal 
surface that separates the attraction field of the ion and that of the target nucleus (see text 
for a more precise definition). According to the results in this paper, the swap-number, 
or number of swaps performed by the electron during the collision appears to be a useful 
parameter in the stopping calculations. In the first place, it turns out that the swap-number 
appropriately labels the electron trajectory, in the sense that it does convey information 
about the extent the electron has interacted with the incoming ion during the collision. At 
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bombarding energies lower than 1 keV the stopping seems to be associated with a fairly 
large number of swaps, whereas, with an increase of the bombarding energy the number of 
swaps relevant to stopping decreases. For bombarding energies larger than 100 keV, however, 
two-swap events appear to be strongly correlated with the stopping process, while both the 
zero and one-swap events become less relevant. Curiously enough, the mean energy- loss 
associated with odd-swap events seems to be a linear function of the bombarding energy, 
whereas those of the even-swaps are not only smaller but also bell-shaped with a maximum 
around 100 keV. 
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